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Abstract 



We present numerical exact results for the ground state and the low-lying excitations 
for the spin-1/2 J1-J2 Heisenberg antiferromagnet on finite square lattices of up to = 40 
sites. Using finite-size extrapolation we determine the ground-state energy, the magnetic 
d • order parameters, the spin gap, the uniform susceptibility, as well as the spin-wave velocity 

^ , and the spin stiffness as functions of the frustration parameter J2 / Ji • In agreement with the 

generally excepted scenario we find semiclassical magnetically ordered phases for J2 < 
Pi ■ and J2 > separated by a gapful quantum paramagnetic phase. We estimate 0.35Ji 

O: and 0.66 Ji. 



PACS codes: 
. 75.10.Jm Quantized spin models 

\ 75.45. +j Macroscopic quantum phenomena in magnetic systems 

l> 

: 1 Introduction 
o . 

^ ■ The spin-1/2 Heisenberg antiferromagnet (HAFM) with nearest-neighbor (NN) Ji and frustrating 
next-nearest-neighbor (NNN) J2 coupling {J1-J2 model) on the square lattice has attracted a great 
deal of interest during the last twenty years (see, e.g., Refs. [1-34] and references therein). The 
^ ' corresponding Hamiltonian reads 

H = Ji^s,-Sj + J2 ^ s, • Sj, (1) 

(iJ) {{i,k)) 

where the sums over {i,j) and {{i, k)) run over all NN and NNN pairs, respectively, counting each 
bond once. In what follows we set Ji = 1. The synthesis of layered magnetic materials [36,37] 
which can be described by the spin-1/2 J1-J2 model has stimulated a renewed interest in this 
model. Another new promising perspective is also opened by the recently discovered layered Fe- 
based superconducting materials [38] which may have a magnetic phase that can be described by 
a J1-J2 model with spin quantum number s > 1/2 [39-41]. 

The J1-J2 model is a canonical model to study quantum phase transitions between semiclas- 
sical magnetically ordered ground state (GS) phases and novel non-magnetic quantum phases. 
Moreover, the J1-J2 model might be also a candidate for a deconfined critical point separating a 
semiclassical magnetic phase from non-magnetic quantum phase [20,30,31,42,43]. 
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Figure 1: The finite square lattice with iV = 40 sites. 



Due to frustration highly efficient quantum Monte-Carlo codes [44], such as the stochastic series 
expansion suffer from the minus sign problem. Therefore, many other approximate methods, e.g., 
the Green's function method [15,22], the series expansion [13,16,18,20,31], the Schwinger boson 
approach [10], the coupled cluster method [12,22,29,30], variational techniques [17,33], the path 
integral quantization [23], the cluster effective- field theory [24, 27], the hierarchical mean-field 
approach [28], the projected entangled pair states method [34] as well as the stochastic state 
selection method [26] were used to find the GS phases of the model. One of the most powerful 
methods is the Lanczos exact diagonalization method which yields numerically exact results for 
finite lattices. It has been succesfully applied to the J1-J2 model in various papers, see, e.g., 
Refs. [2-4,6-9,11,14,21]. In particular, the paper of Schulz and coworkers [4] presenting for the 
first time results for lattices up to iV = 36 sites has set a benchmark for such studies and provides 
detailed information on the model. In addition, these numerical exact results are often used to 
test new approximate methods, see, e.g., Refs. [26,30,35]. 

Due to the progress in the computer hardware and also to the increase in the efficiency in 
programming the Lanczos algorithm, very recently the GS and low-lying excitations of the un- 
frustrated (i.e., J2 = 0) spin-1/2 HAFM [45] and of a spin-1/2 Heisenberg model with ring 
exchange [46] have been calculated for a square lattice with = 40 sites. The largest two- 
dimensional quantum spin model for which the GS has been calculated so far is the spin-1/2 
HAFM on the star lattice with A^ = 42 sites [47]. Note, however, that for higher sectors of 
which are relevant for finite magnetic fields much larger system sizes can be considered, see, e.g., 
Ref. [48]. 

As a result of numerous investigations it seems to be clear now that the J1-J2 model exhibits 
two magnetically long-range ordered phases at small and at large J2 separated by an intermediate 
quantum paramagnetic phase without magnetic long-range order (LRO) in the parameter region 
J^^ < J2 < <^2^ where ~ 0.4Ji and ~ O.6J1. The GS at low J2 < exhibits semiclassical 
Neel magnetic LRO with the magnetic wave vector Qo = (tt, tt). The GS at large J2 > shows 
so-called coUinear magnetic LRO with the magnetic wave vectors Qi = (tt, 0) or Q2 = (0,7r). 
These two collinear states are characterized by a parallel spin orientation of nearest neighbors in 
vertical (horizontal) direction and an antiparallel spin orientation of nearest neighbors in horizontal 
(vertical) direction. 

The nature of the transition between the Neel and the quantum paramagnetic phase as well as 
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the properties of the quantum paramagnetic phase and the precise values of the transition points 
are still under discussion [1-9,12-31,33]. 

In this paper we present new results for the GS and the low-lying excitations of the J1-J2 
model for a finite lattice of = 40 sites shown as in Fig. [H These results are presented in detail, 
on the one hand, as benchmark data for approximate methods. On the other hand, we combine 
our new results with the known results [4] (which we have recalculated) to improve the finite-size 
extrapolation. For all the finite lattices periodic boundary conditions were imposed. 



2 Results of Lanczos exact diagonalization for a finite square 
lattice of = 40 sites 



Although Schulz and coworkers performed the Lanczos diagonalization for the spin-1/2 J1-J2 
HAFM on the finite square lattice of A^ = 36 sites more then 10 years ago the corresponding 
Lanczos diagonalization for the finite lattice of A^ = 40 sites is still a challenging problem. The 
number of basis states on the square lattice in the GS symmetry sector is Uh = 430909650 for 
A^ = 40 compared to Uh = 15804956 for A^ = 36. The GS belongs to a translational symmetry 
with q = (0,0). 

The most important quantities in order to analyze semiclassical GS magnetic ordering are 
appropriate order parameters corresponding to the classical Neel and coUinear LRO. Following 
Schulz et al. [4] we use here the Q-dependent susceptibilities (square of order parameters) defined 
as 



For the magnetic wave vector Qo = (vr, vr) the quantity M^(Qo) is the relevant order parameter 



Qi = (vr, 0) or Q2 = (0, vr) corresponds to the coUinear magnetic LRO present at large J2. 

In tabled] we give the singlet GS energy, Egs{S = 0), the energy of the first triplet excitation 
£'0(5' = 1), the energy of the first singlet excitation Ei{S = 0), as well as the Q-dependent 
susceptibilities M|f(7r, vr), M^(7r, 0). The low- lying energies are also displayed in Fig. [2] where 
more data points are included than shown in table [TJ It is obvious that the first excited singlet 
state becomes very close to the singlet GS near J2 = 0.6. Since both states belong to the same 
lattice symmetry, a level crossing is forbidden due to the von Neumann- Wigner theorem [49] . The 
numerical data shown in the right inset of Fig. [2] clearly give a numerical check of this general 
theorem. The energies of the low-lying states show the typical frustration induced maximum 
known from previous ED calculations for smaller lattices, see, e.g., Refs. [2,4,6]. Since we have 
data for the singlet GS and the first triplet excitation, we can also give results for the spin gap 
At = Eo{S = 1) — Egs{S = 0), see Fig. [31 Note that the spin gap was not presented in the paper 
of Schulz et al. [4]. Hence we compare in Fig. [3] spin gap data for A^ = 32, 36, and 40. It is obvious 
that with increasing J2 the gap starts to grow at J2 ~ 0.4Ji. It reaches a maximum at about 
J2 ~ 0.58 Ji and becomes again small for J2 > 0.6 Ji. This behavior yields an indication of a gapful 
quantum paramagnetic GS phase around J2 = 0.5Ji. However, it is also obvious that there is no 
monotonous finite-size behavior in that parameter region, i.e., a reliable finite-size extrapolation 
of the spin gap around J2 = 0.5Ji is not possible. 

Next we consider the spin-spin correlation functions (sq ■ Sr,). There are altogether 11 different 
correlations functions given in table [2] for the same data points as in table [TJ Graphically the 
variation of some selected correlation functions are shown inf Fig. [H For small J2 the spin-spin 
correlations are quite strong and according to the Neel order we have (sq ■ sr,) < for R connecting 




(2) 




(Qi(2)) with magnetic wave vectors 
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Table 1: Ground state energy Egs{S = 0), first triplet excitation energy Eo{S = 1), first singlet 
excitation energy Ei{S = 0), square of Neel order parameter M^(7r, tt) and square of coUinear 
order parameter M^(7r, 0). 



J2 


Egs{S = 0) 


Eo{S = 1) 




M2,(7r,7r) 




0.0 


-27.09485025 


-26 


83322962 


-26 


3138435 





193923 


0.011956 


0.10 


-25.46460260 


-25 


19683764 


-24 


6669636 





184079 


0.011994 


0.20 


-23.90046918 


-23 


62413648 


-23 


0809149 





171482 


0.012066 


0.30 


-22.42728643 


-22 


13740165 


-21 


5763396 





154648 


0.012212 


0.40 


-21.08836670 


-20 


77332213 


-20 


1877770 





130935 


0.012547 


0.50 


-19.96304839 


-19 


59762345 


-19 


0413364 





097703 


0.013446 


0.55 


-19.51791526 


-19 


11727950 


-18 


8522256 





079043 


0.014457 


0.60 


-19.18368038 


-18 


90110215 


-19 


1762893 





020304 


0.078160 


0.65 


-20.04603255 


-19 


75459787 


-19 


3898477 





013516 


0.093232 


0.70 


-21.05530239 


-20 


75925837 


-20 


2366849 





010120 


0.101552 


0.80 


-23.34020427 


-23 


01691919 


-22 


5277370 





005522 


0.111045 


0.90 


-25.83691287 


-25 


46196339 


-24 


8738325 





002912 


0.115628 


1.00 


-28.43880892 


-28 


00735838 


-27 


2299736 





001700 


0.117879 



sites of sublattices A and B but (sq ■ sr) > for R connecting sites within a sublattice. Increasing 
the frustration leads to a weakening of the spin-spin correlation. This weakening is particularly 
strong for larger separations. For strong frustration around J2 = 0.5Ji all correlation functions 
(sq -Sr.) except for R = (1, 0); (1,1); (2, 0) are very small which is an indication for a magnetically 
disordered phase. Beyond J2 ~ O.6J1 the increasing strength of correlation functions for lattice 
vectors R connecting sites within the same sublattice A (or B) indicates the emerging collinear 
LRO. The corresponding order parameters defined in Eq. ([2]) and listed in table [1] are shown 
in Fig. [5l The suppression of magnetic order around J2 = 0.5Ji is again obvious. Although 
there is no level crossing (see the discussion above) we observe in Figs. H] and [5] a sharp change 
of the magnetic quantities near J2 = O.6J1, i.e., at that point where both singlet levels become 
very close to each other (see Fig. [21 right inset). This behavior might be a hint on a first-order 
transition between the phase with collinear LRO and the magnetically disordered phase, see also, 
e.g., Refs. [4,13,16,22,30]. 

Another interesting point is the breakdown of the Marshall-Peierls sign rule (MPSR) at J2 ~ 
0.4Ji. Writing the GS as = ^„c„|?t,), where \n) is an Ising basis state of typical form 
I TTitii ■■•), the MPSR determines the sign of the coefficients [50]. The MPSR has been 
proved exactly for bipartite lattices and arbitrary site spins by Lieb, Schultz and Mattis [51]. 
As pointed out in several papers the knowledge of the sign of the c„ is of great importance in 
different numerical methods, e.g. for the construction of variational wave functions [7,17,33,52], 
in quantum Monte-Carlo methods (which suffer from the sign problem in frustrated systems [53]) 
and also in the density matrix renormalization group method, where the application of the MPSR 
has substantially improved the method in a frustrated spin system [54]. 

For the J1-J2 model on the square lattice the violation of the MPSR was considered as one 
signal of destruction of magnetic LRO [7,8,11,12,17,33]. It was found that the MPSR (proved 
exactly only for J2 = 0) survives till about J2 ~ 0.3 . . . 0.4Ji. The exact Lanczos results for the 
ground state presented in Refs. [7,8] and [11] were, however, restricted to systems of up to iV = 24 
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Table 2: Spin-spin correlation function (sq ■ Sr) (respectively, (sq ■ Sj), see Fig. [T]) 





(so ■ Si) 




(so ■ S3) 


(So ■ S4) 


(so ■ S5) 


(so ■ se) 


J2 


R= (1,0) 


R= (2,0) 


R 


= (3,0) 


R=(3,l) 


R= (2,1) 


R=(l,l) 


0.0 


-0.338685 


0.185133 


-0 


166989 


0.164100 


-0.179037 


0.207198 


0.10 


-0.338313 


0.175458 


-0 


154968 


0.151833 


-0.166128 


0.200046 


0.20 


-0.336852 


0.163089 


-0 


140055 


0.136281 


-0.149688 


0.190488 


0.30 


-0.333423 


0.146655 


-0 


120879 


0.115776 


-0.127803 


0.176943 


0.40 


-0.326073 


0.123909 


-0 


095208 


0.087441 


-0.097056 


0.156171 


0.50 


-0.310788 


0.093354 


-0 


061983 


0.048666 


-0.053772 


0.122499 


0.55 


-0.298596 


0.077397 


-0 


045300 


0.027102 


-0.028947 


0.099315 


0.60 


-0.129468 


0.145560 


-0 


032355 


-0.127482 


0.034386 


-0.183879 


0.65 


-0.096630 


0.168312 


-0 


027747 


-0.156408 


0.033291 


-0.236841 


0.70 


-0.077052 


0.181365 


-0 


024825 


-0.170640 


0.030228 


-0.265914 


0.80 


-0.050388 


0.195918 


-0 


020247 


-0.186357 


0.025323 


-0.301704 


0.90 


-0.034740 


0.202785 


-0 


016548 


-0.193845 


0.021465 


-0.320247 


1.00 


-0.026295 


0.206295 


-0 


013959 


-0.197463 


0.018474 


-0.329190 





(So ■ Sg) 


(so • Sio) 


(so ■ Sii) 


(so ■ Sie) 


(Sq ■ S24) 


J2 


R=(l,2) 


R=(l,3) 


R= (2,2) 


R= (2,3) 


R = (4, -2) 


0.0 


-0.177165 


0.157032 


0.161478 


-0.156222 


0.150825 


0.10 


-0.164196 


0.144276 


0.148815 


-0.142812 


0.137613 


0.20 


-0.147774 


0.128166 


0.132678 


-0.125913 


0.120888 


0.30 


-0.126123 


0.107007 


0.111225 


-0.103674 


0.098736 


0.40 


-0.096186 


0.077940 


0.081246 


-0.072885 


0.067683 


0.50 


-0.055281 


0.039045 


0.040050 


-0.031119 


0.024636 


0.55 


-0.032856 


0.018723 


0.018048 


-0.009279 


0.001776 


0.60 


0.015885 


-0.109425 


0.107499 


0.004956 


0.092322 


0.65 


0.017844 


-0.138945 


0.137898 


0.002280 


0.124962 


0.70 


0.017595 


-0.155349 


0.155199 


0.000924 


0.142587 


0.80 


0.016110 


-0.173610 


0.174066 


0.000216 


0.159999 


0.90 


0.014244 


-0.182238 


0.182661 


0.000291 


0.166980 


1.00 


0.012522 


-0.186480 


0.186960 


0.000333 


0.170301 
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-0.45 



-0.5 
-0.55 
-0.6 

z 

LLJ 

-0.65 



-0.7 
-0.75 



-0.8 

0.2 0.4 0.6 0.8 1 

J2/J1 

Figure 2: Singlet GS energy, Ecs{S = 0), the energy of the first triplet excitation, Eo{S = 1), and 
the energy of the first singlet excitation, Ei{S = 0), versus J2 for the square lattice with = 40 
sites. 





J2/J1 



Figure 3: Spin gap = Egs{S = 1) - ^o('^ = 0) versus J2 for N = 32, 36, and 40. 
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Figure 5: Square of order parameters M^(Q) for Q = (tt, tt) and Q = (0, tt), see Eq. ([2]). 
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J2/J1 



Figure 6: Validity of the MPSR: The curve shows the weight of the states fulfiUing the MPSR, 
i.e., |c„|^ (where the sum ^ runs over those states only, for which the sign of c„ fulfills the 
MPSR). The arrows indicate those J2 values for TV = 40, where the weight of the states which do 
not fulfill the MPSR is 10"^^ lO'^ and 10"^ respectively. 

sites. Here we present data for = 32, 36, and 40 in Fig. [61 From Fig. [6] it is obvious that 
the MPSR is valid in very good approximation even till J2 ~ 0.4Ji. The MPSR starts to be 
significantly violated beyond J2 = 0.4Ji, where the Neel LRO breaks down. 



3 Finite-size extrapolation 

We use now our new data for = 40 together with the known data for N = 20, 32, 36 to 
perform a finite-size extrapolation. Note that we do not include the data for the lattice with 

= 16 sites. This lattice has an extra symmetry because it is equivalent to a hypercube in 4 
dimensions. It was argued that, therefore, the N = 16 lattice exhibits an anomalous behavior [4]. 
The finite-size extrapolation rules for the two-dimensional HAFM are well known [55-59]. The 
extrapolation of the magnetic order parameters yields an estimate of the transition points between 
the semiclassically ordered phases (Neel and coUinear) and the magnetically disordered quantum 
phase. In addition, since the scaling behavior is related to the low-energy degrees of freedom of 
the model, via the extrapolation procedure one can extract the J2 dependence of other quantities 
such as the spin-wave velocity c and the spin stiffness p. It has been demonstrated that Lanczos 
data for finite lattices of up to N=36 sites can provide quite accurate data for these quantities in 
the case of unfrustrated lattices [4,58,60]. 

In the seminal paper of Schulz et al. [4] the finite-size extrapolation procedure as well as some 
particular problems of the extrapolation appearing for the J1-J2 model were discussed in great 
detail. We follow here the lines of that paper, but do not repeat a detailed discussion of the 
extrapolation scheme. 

We define the staggered magnetization, i.e., the Neel order parameter, as mo = 2 limTv^oo Mjsflir, vr). 
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Figure 8: Ground-state energy per site. Main panel: Extrapolated value cq (see Eq. ([5]). Inset: 
Data for = 20, 32, 36, and = 40 around the maximum. 
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The finite-size behavior of M7v(vr, tt) is given by [4,55-57] 

,.2. X 1 2A 0.62075c \ 

MUn,n) = -ml^l + —-^ + ...J (3) 

where c is the spin-wave velocity and p is the spin stiffness. The corresponding order parameter 
for the colhnear LRO at large J2 is defined as [4] nii = -\/8 liniTv^oo AfAr(7r, 0). The finite-size 
behavior of Mn^h, 0) is 

, ,9 / N 1 9 const. , 
M2(7r,0) = -m2 + — ^ + ... . (4) 



The results for mo and mi are shown in Fig. [71 As expected we find a magnetically disordered 
GS phase around J2 = 0.5Ji. The transition points are determined to = 0.35Ji and = 
O.66J1, i.e., the range of the magnetically disordered GS phase obtained here is slightly larger 
than predicted, e.g., from series expansion [13,16,18,20] or coupled cluster approach [30]. For the 
unfrustrated square lattice HAFM ( J2 = 0) we obtain mo = 0.621 which is in good agreement with 
corresponding results obtained by other methods, e.g., 3rd order spin-wave theory (mo = 0.6138) 
[63], quantum Monte Carlo (mo = 0.6140) [62] or coupled cluster method (mo = 0.6205) [61]. 
The finite-size behavior of the GS energy is given by 

^^ = eo- 1.4372^ + --- . (5) 

The extrapolated energy cq is shown in Fig. [HI Because of an irregular finite-size behavior near 
the maximum in cq (see inset in Fig. [Hj) the extrapolation becomes unreliable around J2 = 0.5Ji 
(this parameter region is indicated by the dotted line in the main panel of Fig. [8]). Nevertheless, 
one can speculate that the kink in the extrapolated GS energy near J2 = O.6J1 might be another 
hint for a first-order transition at between the semiclassical collinear phase and the quantum 
paramagnetic phase. The extrapolated GS energy, co = —0.6701, for J2 = is again in very good 
agreement with corresponding results obtained by other methods, e.g., 3rd order spin- wave theory 
(eo = —0.66931) [63], quantum Monte Carlo (co = —0.66944) [62] or coupled cluster method 
(eo = -0.66936) [61]. 

We can also perform a finite-size extrapolation of the gap to the first triplet excitation (spin 
gap) Aj'(A^) = Eo{S = 1) — Egs{S = 0). The corresponding formula is 

At{N) = At + ^ + --- . (6) 

The results for are shown in Fig. [HI In the magnetically ordered Neel and collinear phases 
gapless Goldstone modes exist and consequently the spin gap should vanish. The finite values of 
the extrapolated gap for J2 < J2^ and J2 > ^2^ give a hint of the limits of precision of the finite- 
size extrapolation of the spin gap. Nevertheless, it becomes obvious that there is a significant 
increase of At in the quantum paramagnetic phase < J2 < J2^ ■ Note that for larger J2 the 
energy scale of the system becomes proportional to J2 that explains the increase in with J2 
for J2 > 0.8. A finite spin gap for < J2 < J2^ would be consistent with the findings of many 
papers which suggest that the quantum paramagnetic phase exhibits valence-bond order, see, e.g., 
Refs. [20] and [30] and references therein. Note, however, that our finding is in contrast to recent 
DMRG studies on frustrated odd-leg ladders [25] where arguments for a gapless spectrum of the 
J1-J2 model on the square lattice have been given. 

Next we consider the uniform susceptibility x which is related to the spin gap by x{^) = 
NAt{N). For the finite-size extrapolation of the uniform susceptibility we use [4] 

x{N) = x + ^ + --- ; xiN)=NAT{N). (7) 
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Figure 9: Spin gap At for = 40 and extrapolated to ^ oo (see Eq. (|6])). 
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J2/J1 



Figure 11: Spin wave velocity c and spin stiffness p. For better comparison p is multiplied by 10. 

The results for xi^ = 40) and the extrapolated value x are shown in Fig.[Tni In a gapful quantum 
paramagnetic phase, i.e., for < J2 < ^2^, x should vanish. Indeed, we observe a vanishing of x 
around J2 = 0.5Ji. However, the J2 parameter region where x is zero is clearly smaller than that 
where the order parameters mo and mi vanish. At J2 = we obtain x = 0.0674 which is in good 
agreement with data obtained by other methods, e.g., quantum Monte Carlo (x = 0.0669) [62], 
series expansion (x = 0.0659) [64], third-order spin-wave theory (x = 0.06291) [63], or coupled 
cluster method (x = 0.065) [65]. 

So far we have considered quantities which can be calculated directly for each finite lattice and 
for which the corresponding value for — ^ 00 is the leading term in a finite-size extrapolation 
formula, see Eqs. ([3]), (j4]), ([5]), ([6]), and ([7j). A test of the accuracy the finite-size extrapolation is 
given by the comparison with best available data for the unfrustrated limit (J2 = 0), see above. In 
addition to the direct calculation of certain magnetic quantities the finite size-extrapolation allows 
also an indirect determination of the spin-wave velocity c and the spin stiffness p which enter the 
finite size-extrapolation formulas ([3]) and ([5]) as prefactors of the leading finite-size corrections. 
First we can determine the spin- wave velocity c via Eq. ([5]). Using this result and also m^ we 
find then the stiffness pi oc ttIqC via Eq. Q. Since pi is proportional to ttIq, it vanishes at the 
same point J2^. Alternatively, we can use the hydrodynamic relation x = p/c^, see, e.g., Refs. [66] 
and [62], to determine the stiffness p2 = xc^- We show the spin- wave velocity and both values of 
the stiffness in Fig. [TTl Although, pi obtained via Eq. ([3]) and p2 obtained via the hydrodynamic 
relation are quite different the qualitative behavior of both is similar to that obtained by direct 
calculation of p, e.g., using Schwinger boson [10] or coupled cluster approach [30]. The spin-wave 
velocity c also decreases with growing J2. However, c remains finite at the magnetic- nonmagnetic 
transition at as predicted from the non-linear sigma model [66]. At J2 = we get c = 1.287, 
pi = 0.1246 and p2 = 0.1115. These values are significantly lower than corresponding values 
obtained by quantum Monte Carlo method [62], 3rd order spin- wave theory [63] or coupled cluster 
method [61,67]. Hence we may argue that the indirect determination of magnetic quantities via 
the prefactors of the leading finite-size corrections is less accurate. 
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4 Summary 



In this paper we have presented detailed information on the GS and the low-lying excitations for 
the spin-1/2 J1-J2 HAFM on a finite square lattice of = 40 sites. These data for the GS energy, 
the lowest triplet and singlet excitations, the spin-spin correlation functions and the magnetic 
order parameters may serve as benchmarks for approximate methods. 

Including data for finite lattices of = 20, 32, 36 sites we have performed finite-size extrap- 
olations of several magnetic quantities. Using the well-known extrapolation formulas for the GS 
energy, the order parameters, the spin gap, and the uniform susceptibility we have determined 
these quantities for N 00. To estimate the accuracy of the extrapolated quantities we have 
compared them with best available results for the unfrustrated limit J2 = 0, where the high-order 
spin-wave theory and the quantum Monte Carlo method work well. At J2 = the extrapolated 
GS energy is in excellent agreement with spin-wave theory and quantum Monte Carlo results. The 
deviation of the Neel order parameter obtained by the finite-size extrapolation from the spin-wave 
and quantum Monte Carlo results is only about 1%. 

Besides these quantities which have been determined directly for each finite lattice, in addition, 
we have used the formulas of finite-size extrapolation to determine the spin-wave velocity as well 
as the spin stiffness which appear in the prefactors of the leading finite-size corrections of the GS 
energy and the Neel order parameter. This 'indirect way' to determine spin-wave velocity and 
stiffness, however, yields only qualitative agreement with known results. 

From the extrapolated magnetic order parameters we find the transition point from the mag- 
netically ordered Neel phase to the gapful quantum paramagnetic phase to ~ 0.35Ji and 
transition point between the magnetically ordered coUinear phase and the quantum paramagnetic 
phase to J 2^ O.66J1. 
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